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METHODS AND SYSTEMS FOR DETECTING 
COMPONENTS OF PLAQUE 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the priority of U.S. Provisional Patent 
Application Serial Number 60/398,341 filed July 23, 2002, which is hereby incorporated by 
referenced in its entirety. 

BACKGROUND OF THE INVENTION 

[0002] This invention relates to imaging systems, and more particularly to 
methods and systems for detecting components of plaque. 

[0003] Recent advancements in computed tomography (CT) technology have 
resulted in faster scanning speeds, larger coverage with multiple detector rows, and thinner 
slices. However, energy resolution in measurements obtained with these systems is still a 
missing piece. The wide x-ray photon energy spectrum output from the x-ray source and the 
lack of energy resolution from CT detection systems preclude energy discrimination CT. 

[0004] The X-ray attenuation through a given object is strongly dependent 
on the incident x-ray photon energy. This physical phenomenon manifests itself in an image 
as beam-hardening artifacts, such as, non-uniformity, shading, and streaks. Some beam- 
hardening artifacts can be easily corrected, but other beam-hardening artifacts may be more 
difficult to correct. In general, known methods to correct beam hardening artifacts include 
water calibration, which includes calibrating each CT machine to remove beam hardening 
from materials similar to water, and iterative bone correction, wherein bones are separated in 
a first-pass image reconstruction then the measurements are corrected for the beam hardening 
resulting from the bones in a second-pass reconstruction. However, beam hardening from 
materials other than water and bone, such as metals and contrast agents, may be difficult to 
correct without added processing of the measured data. In addition, even with the above 
described correction methods, conventional CT does not provide quantitative image values. 
Rather, the same material at different locations within the image often shows a different CT 
number. 
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[0005] Another drawback of conventional CT is the reduced level of 
material characterization. For example, a highly attenuating material with a low density can 
result in the same CT number in the image as a less attenuating material with a high density. 
Thus, there is little or no information about the material composition of a scanned object 
based solely on the CT number. 

BRIEF DESCRIPTION OF THE INVENTION 

[0006] In one aspect, a method is provided. The method includes detecting 
components of plaque using a multi-energy computed tomography (MECT) system. 

[0007] In another aspect, a method for detecting components of plaque is 
provided. The method includes generating information regarding projection data of phantoms 
by using an MECT system, and obtaining the components of the plaque from the information. 

[0008] In yet another aspect, an MECT system is provided. The MECT 
system includes at least one radiation source configured to transmit x-rays that intersect an 
object; at least one detector configured to detect the x-rays; a controller coupled to the 
detector; and a computer configured to instruct the MECT system to detect components of 
plaque. 

[0009] In still another aspect, a computer readable medium encoded with a 
program is provided. The program is configured to instruct a computer to detect components 
of plaque within an object that is scanned using a multi-energy tomography (MECT) system. 

[0010] In still another aspect, a computer encoded with a program is 
provided. The program is configured to instruct an MECT system to detect components of 
plaque within an object that is scanned using the MECT system. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] Figure 1 is a pictorial view of a multi-energy computed tomography 
(MECT) imaging system in which a method for detecting components of plaque is 
implemented. 
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[0012] Figure 2 is a block schematic diagram of the MECT system 
illustrated in Figure 1 . 

[0013] Figures 3, 4, 5, and 6 illustrate a flowchart of an embodiment of the 
method for detecting components of plaque. 

[0014] Figure 7 illustrates a plurality of images. 

[00 1 5] Figure 8 illustrates another plurality of images. 

DETAILED DESCRIPTION OF THE INVENTION 

[0016] It is believed that progression of arteriosclerosis causes a build-up of 
fat-laden deposits on vessel walls, causing progressive narrowing in the vessel, such as in a 
carotid or a coronary artery. Eventually, lumen or blood flow within the vessel is reduced to 
such a level that a tissue, such as, a heart muscle or brain tissue, is starved of oxygen-carrying 
blood and the tissue dies, resulting in a heart attack or stroke. Comparatively, in 
atherosclerosis, a majority of patients with acute cardiac disease exhibit an inflammatory 
process within the vessel wall itself. Low-density lipoproteins (LDLs) accumulate in the 
vessel wall and attract immune system cells into the vessel wall as well. The immune system 
cells ingest the modified LDLs, giving rise to fatty droplets, which constitute a lipid core of a 
plaque. The inflammatory cells promote growth of the plaque and development of a fibrous 
cap composed of collagen fibers formed by smooth muscle cells over the lipid core. When 
the inflammatory cells flare, immune system cells secrete enzymes that degrade collagen and 
prevent the development of new collagen fibers to repair cap damage. The weakening of the 
cap may continue until a fissure occurs and the blood of the lumen intermingles with the lipid 
core, rich in proteins that foster blood coagulation. As a result, a clot forms and the vessel 
may be occluded, resulting in an ischemic event. If the clot is not occlusive, scar tissue 
develops and the plaque is enlarged, giving rise to a chronic condition. 

[0017] In coronary atherosclerosis, thrombi, as well as calcified, 
intermediate, and soft plaques could be present. Plaque composition is believed to be 
indicative of the risk of acute coronary syndromes. Soft plaque includes a high lipid 
concentration, a fibrous thin cap and inflammatory cells. Intermediate plaques include 
fibrous tissue, lower concentrations of lipids and inflammatory cells. Calcified plaques 
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include a high concentration of calcium. If soft plaque is present, there may be an increased 
risk of a stroke or a heart attack. It is believed that thrombi have a density less than 20 
Hounsfield units (HU), soft plaques have a density less than 50 HU, intermediate plaques 
have a density between 50 and 120 HU, and calcified plaques have a density over 120 HU. 
Hounsfield units, also known as CT numbers, are used to control the brightness of a 
corresponding pixel on a cathode ray tube display and reflect the consistency of the different 
types of plaques. However, it is challenging to identify and discriminate plaques with low or 
no concentration of calcium. 

[0018] In some CT imaging system configurations, an x-ray source projects a 
fan-shaped beam which is collimated to lie within an X-Y plane of a Cartesian coordinate 
system and generally referred to as an "imaging plane". The x-ray beam passes through an 
object being imaged, such as a patient. The beam, after being attenuated by the object, 
impinges upon an array of radiation detectors. The intensity of the attenuated radiation beam 
received at the detector array is dependent upon the attenuation of an x-ray beam by the 
object. Each detector element of the array produces a separate electrical signal that is a 
measurement of the beam attenuation at the detector location. The attenuation measurements 
from all the detectors are acquired separately to produce a transmission profile. 

[0019] In third generation CT systems, the x-ray source and the detector 
array are rotated with a gantry within the imaging plane and around the object to be imaged 
such that the angle at which the x-ray beam intersects the object constantly changes. A group 
of processed x-ray attenuation measurements, i.e. projection data, from the detector array at 
one gantry angle is referred to as a "view". A "scan" of the object comprises a set of views 
made at different gantry angles, or view angles, during one revolution of the x-ray source and 
detector. 

[0020] In an axial scan, the projection data is processed to construct an 
image that corresponds to a two dimensional (2D) slice taken through the object. One 
method for reconstructing an image from a set of projection data is referred to in the art as the 
filtered back projection technique. This process converts the attenuation measurements from 
a scan into integers called Hounsfield units. 
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[0021] To reduce the total scan time, a "helical" scan may be performed. To 
perform a "helical" scan, the object is moved into the scanner while the data for the 
prescribed number of slices is acquired. Such a system generates a single helix from an x-ray 
fan beam helical scan. The helix mapped out by the x-ray fan beam yields projection data 
from which images in each prescribed slice may be reconstructed. 

[0022] Reconstruction algorithms for helical scanning typically use helical 
weighing algorithms that weight the collected data as a function of view angle and detector 
channel index. Specifically, prior to a filtered backprojection process, the data is weighted 
according to a helical weighing factor, which is a function of both the gantry angle and 
detector angle. The helical weighting algorithms also scale the data according to a scaling 
factor, which is a function of the distance between the x-ray source and the object. The 
weighted and scaled data is then processed to generate CT numbers and to construct an image 
that corresponds to a 2D slice taken through the object. 

[0023] As used herein, an element or step recited in the singular and 
preceded with the word "a" or "an" should be understood as not excluding plural said 
elements or steps, unless such exclusion is explicitly recited. Furthermore, references to "one 
embodiment" of the present invention are not intended to be interpreted as excluding the 
existence of additional embodiments that also incorporate the recited features. 

[0024] Also as used herein, the phrase "reconstructing an image" is not 
intended to exclude embodiments of the present invention in which data representing an 
image is generated but a viewable image is not. However, many embodiments generate (or 
are configured to generate) at least one viewable image. 

[0025] Referring to Figures 1 and 2, a Multi-Energy multi-slice scanning 
imaging system, for example, a Multi-Energy computed tomography (MECT) imaging system 
10, is shown as including a gantry 12 representative of a "third generation" CT imaging 
system. Methods for detecting components of plaque are implemented using MECT system 
10. Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a detector 
array 18 on the opposite side of gantry 12. Detector array 18 is formed by a plurality of 
detector rows (not shown) including a plurality of detector elements 20 which together sense 
the projected x-rays that pass through an object, such as a medical patient 22. Each detector 
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element 20 produces an electrical signal that represents the intensity of an impinging x-ray 
beam and hence can be used to estimate the attenuation of the beam as it passes through 
patient 22. During a scan to acquire x-ray projection data, the components mounted on 
circular portion of gantry 12 therein rotate about a center of rotation 24. Figure 2 shows only 
a single row of detector elements 20 (i.e., a detector row). However, multi-slice detector 
array 18 includes a plurality of parallel detector rows of detector elements 20 such that 
projection data corresponding to a plurality of quasi-parallel or parallel slices can be acquired 
simultaneously during a scan. Moreover, detector array 18 may be an area detector 
facilitating large imaging coverage on the patient 22 for each rotation of the circular portion 
of gantry 12. 

[0026] Rotation of components in gantry 12 and the operation of x-ray 
source 14 are governed by a control mechanism 26 of MECT system 10. Control mechanism 
26 includes an x-ray controller 28 that provides power and timing signals to x-ray source 14 
and a gantry motor controller 30 that controls the rotational speed and position of components 
in gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples analog 
data from detector elements 20 and converts the data to digital signals for subsequent 
processing. An image reconstructor 34 receives sampled and digitized x-ray data from DAS 
32 and performs high-speed image reconstruction. The reconstructed image is applied as an 
input to a computer 36 which stores the image in a storage device 38. Image reconstructor 
34 can be specialized hardware or computer programs executing on computer 36. 

[0027] Computer 36 also receives commands and scanning parameters from 
an operator via console 40 that has a keyboard. An associated cathode ray tube display 42 
allows the operator to observe the reconstructed image and other data from computer 36. The 
operator supplied commands and parameters are used by computer 36 to provide control 
signals and information to DAS 32, x-ray controller 28, and gantry motor controller 30. In 
addition, computer 36 operates a table motor controller 44, which controls a motorized table 
46 to position patient 22 in gantry 12. Particularly, table 46 moves portions of patient 22 
through gantry opening 48. 

[0028] In one embodiment, computer 36 includes a device 50, for example, a 
floppy disk drive or CD-ROM drive, for reading instructions and/or data from a computer- 
readable medium 52, such as a floppy disk or CD-ROM. In another embodiment, computer 

-6- 



127005 

Express Mail No.: EV299866125US 



36 executes instructions stored in firmware (not shown). Computer 36 is programmed to 
perform functions described herein, and as used herein, the term computer is not limited to 
just those integrated circuits referred to in the art as computers, but broadly refers to 
controllers, processors, microcontrollers, microcomputers, programmable logic controllers, 
application specific integrated circuits, and other programmable circuits, and these terms are 
used interchangeably herein. MECT system 10 is an energy-discriminating computed 
tomography (EDCT) system in that MECT system 10 is configured to be responsive to 
different incident x-ray spectra. This can be accomplished with a conventional third 
generation CT system by acquiring projection data sequentially using different x-ray tube 
voltages. For example, two scans are acquired either back to back or interleaved in which the 
tube operates at 80kVp and 160kVp potentials, for example, generating a low and high energy 
spectrum respectively. Alternatively, special filters are placed between x-ray source 14 and 
detector array 18 such that different detector rows collect projections of different x-ray energy 
spectrum. Yet another embodiment is to use energy sensitive detectors such that each x-ray 
photon reaching the detector is recorded with its photon energy. Although the specific 
embodiment mentioned above refers to a third generation CT system, the methods described 
herein equally apply to fourth generation CT systems (stationary detector - rotating x-ray 
source) and fifth generation CT systems (stationary detector and x-ray source), as well as 
other higher generation CT systems. 

[0029] EDCT can lessen or eliminate the problems, such as lack of energy 
discrimination or material characterization, associated with some CT systems altogether. In 
the absence of object scatter, MECT system 10 may be used to separately detect two regions 
of the incident photon energy spectrum, the low-energy and the high-energy portions of the 
incident x-ray spectrum. The behavior at any other energy can be derived based on the signal 
from the two energy regions. This phenomenon is driven by the fundamental fact that in the 
energy region where medical CT is interested, two physical processes dominate the x-ray 
attenuation (1) Compton scatter and (2) the Photoelectric effect. In order to characterize the 
behavior of an object causing attenuation of the x-ray beam, two independent parameters are 
measured. Thus, detected signals from the two energy regions provide sufficient information 
to resolve the energy dependence of the object being imaged; hence, the composition of the 
material can be characterized. 
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[0030] The data analysis used in EDCT includes Compton and photoelectric 
decomposition and/or Basis material decomposition (BMD). In Compton and photoelectric 
decomposition, a pair of images is generated with MECT 10, which separately present the 
attenuation from the Compton and photoelectric processes - instead of obtaining one image 
characterizing the overall attenuation coefficient in the reconstructed CT image. Also, a slight 
modification in the processing allows the generation of images representing density and 
effective atomic number. The BMD method is based on the concept that the x-ray attenuation 
of any given material in the energy range for MECT can be represented by assuming a density 
mixture of other two known materials. These two materials are called the Basis Materials. 
Using BMD, two reconstructed images are obtained, each image representing the equivalent 
density of one of the basis materials. Since density is independent of x-ray photon energy, 
these images are relatively free of beam hardening artifacts. Additionally, the basis material is 
chosen to target a material of interest, thus enhancing the image contrast. 

[003 1 ] It should be noted that in order to optimize a multi-energy CT system 
not implementing energy discrimination with photon counting, the larger the energy 
separation in the x-ray spectra, the better the image quality. Also, the photon statistics in these 
two energy regions should be comparable, otherwise the energy region with reduced 
statistical information will dominate the noise in the reconstructed image. 

[0032] There are different methods to obtain dual energy measurements: (1) 
scan with two distinctive energy spectra, (2) detect photon energy according to penetration 
depth at the detector, or (3) photon-counting with energy discrimination. Photon counting 
provides clean spectra separation and an adjustable energy separation threshold for balancing 
photon statistics. 

[0033] Figures 3, 4, 5, and 6 illustrate a flowchart of an embodiment of a 

method for detecting components of a plaque, such as atherosclerotic plaque. The method is 

executed by computer 36 upon receiving images reconstructed by image reconstructor 34. 

The method includes detecting 66 components of a plaque, such as calcium, lipid, and smooth 

muscle cells, developed in patient 22 and includes quantification 68 of the components of 

plaque. MECT system 10 is used to detect the components. The components are detected by 

generating information regarding densities of constituent components of the plaque and 

characterizing the components of the plaque from the information. The information may be 
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generated by simulating the data collection of projection data obtained with an object, such 
as, for instance, a phantom, to allow generation of reconstructed images of the phantom by 
MECT system 10. Material decomposition basis materials, such as water and iodine, of the 
phantom are assumed. The phantom includes regions with known material properties and 
known densities of water and iodine. For instance, the phantom includes a region with an 
iodine solution of known concentration, and thus known water and iodine densities. As 
another example, the phantom includes a region with an iodine solution of a different known 
concentration. 

[0034] To generate the information, a first and the second sets of projection 
data are generated by spanning the entire space of iodine and water concentrations of interest 
within the phantom. The first set of projection data is generated by performing 70 a 
simulation of placing the phantom in a scanning field of MECT system 1 0 and scanning 72 
the phantom at a first energy spectrum using MECT system 10. As an example, the first 
energy spectrum is a high or a low energy spectrum. Therefore, for known water and iodine 
densities, for a known material, projection data that correspond to the first energy spectrum 
are obtained. The second set of projection data of the phantom is obtained 74 in a similar 
manner as obtaining the first set of projection data except that the second set of projection 
data is generated by simulating the scanning of the phantom with a second energy spectrum. 
The second energy spectrum is a low energy spectrum if the first energy spectrum is a high- 
energy spectrum. Conversely, the second energy spectrum is a high-energy spectrum if the 
first energy spectrum is a low-energy spectrum. Once the simulated data acquisition is 
completed for the first and second energy spectrum scans, for all the known different iodine 
concentrations, a specific material basis decomposition look-up table is formed 76. If iodine 
concentrations are used, the look-up table maps known projection data generated from 
measurements of water and iodine in the phantom to allow generation of density distributions 
in first and second sets of images, where the first set of images is generated from the first set 
of projection data and the second set of images is generated from the second set of projection 
data. After the construction of the look-up table, the components of the plaque are obtained 
by scanning 78 an object, such as patient 22, of unknown material to generate projection data 
and using 80 the look-up table to allow generation of images of the density of the basis 
materials. As an example, the look-up table is reversely mapped to obtain material 

distributions within patient 22 from densities of material distributions within the phantom. 
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[0035] A detailed description of the step of using the look-up table to allow 
generation of images of the density of the basis materials has been described in an application 
titled "Method and Apparatus For Generating a Density Map Using Dual-Energy CT", having 
attorney docket number 127076-2, and filed on July 22, 2003, which is incorporated by 
reference herein in its entirety. The detailed description of the step of using the look-up table 
to allow generation of images of the density of the basis materials has also been described in a 
provisional application having serial number 60/397,658 and filed on July, 23, 2002, which is 
also incorporated by reference herein in its entirety. 

[0036] The components of the plaque are obtained by placing patient 22 in 
the scanning field of MECT system 10 and scanning patient 22 at the first and second energy 
spectrums using MECT system 10 to obtain reconstructions of the density of the basis 
materials within patient 22. In an alternative embodiment, the method includes scanning a 
first small cross section of patient 22. The first small cross section includes cells or elements 
that exist in a lipid core of the plaque but do not exist in the smooth muscles of the plaque. In 
yet another alternative embodiment, the method includes scanning a second small cross 
section of patient 22. The second small cross section includes the plaque. It is noted that 
additional energy levels can be used to scan the phantom or the object to generate the look-up 
table. 

[0037] The method further includes quantifying 68 the components of the 
plaque. The method quantifies the components of the plaque by calculating composition 
density distributions, such as lipid burden or calcium burden, of the plaque and calculating 
total plaque burden. As an example, the total plaque burden is calculated by a thresholding 
method. The thresholding method distinguishes pixels that have a density greater than a 
number y. The pixels are pixels of images obtained by using MECT system 10. The total 
plaque burden in an organ of patient 22 is a weighted sum of densities greater than y of pixels 
of an image of the organ. The composition density distributions are calculated by a method 
similar to that used to calculate the total plaque burden. For instance, the lipid burden is 
calculated by distinguishing pixels of images obtained using MECT system 10 that have a 
density greater than a number v. The total lipid burden in an organ of patient 22 is a weighted 
sum of densities greater than v of pixels of an image of the organ. In an alternative 
embodiment, the quantities of the components are combined with measurements of geometry 
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of lumen of patient 22 to assess degree and severity of stenosis. Other techniques for 
quantifying components of plaque are described in "Coronary Artery Calcium: Accuracy and 
Reproducibility of Measurements with Multi-Detector Row CT-Assessment of Effects of 
Different Thresholds and Quantification Methods 1," Cheng Hong, Kyongtae T. Bae and 
Thomas K. Pilgram, Radiology 2003, 227: 795-801; "Noninvasive quantification of coronary 
artery calcification: Methods and prognostic value/' Sandra S. Halliburton, Arthur E. 
Stillman, Richard D. White, Cleveland Clinic Journal Of Medicine Volume 69, Supplement; 
and in "CT and MR for coronary Angiography & plaque imaging," Fayad et aL, Circulation, 
Oct. 2002, pp. 2026-2033. 

[0038] The method also includes instructing 104 MECT system 10 to 
perform additional scans of patient 22 at different times and repeating step 66 shown in 
Figure 3 each time a scan is performed. The method includes waiting 106 for a user to 
administer a contrast agent, such as non-ionic Iodine, Gadolinium, a blood pool contrast 
agent, and/or a contrast agent that adheres to specific molecules. The contrast agent is 
administered to accentuate the contrast in at least one of lipid-avid agents, plaque specific 
antigens, and/or plaque cells of patient 22. The method further includes repeating 107 step 
66. 

[0039] The method also includes waiting 108 for the user to administer a 
temperature-sensitive contrast agent to accentuate an inflamed plaque. Examples of the 
temperature-sensitive contrast agents include agents containing a lanthanide metal, such as 
Ce, Pr, Nd, Sm, Eu, Gd, Db, Dy, Ho, Er, Tm, and Yb. The method includes repeating 109 
step 66. 

[0040] As identified in Figure 5, the method further includes displaying 110 

at least one of a 2-dimensional (2D) and a 3-dimensional (3D) image of the components of 

the plaque on a wall of an organ, such as an intestine or blood vessel of patient 22, and 

enabling viewing, from a viewpoint, a volume of the plaque in the 3D image. The method 

performs volume rendering for viewing the 3D image or using visualization methods to 

traverse through the 3D image from at least one viewpoint. Volume rendering is a technique 

used for visualizing sampled functions in three spatial dimensions. Volume rendering is 

applied to medical imaging where volume data is available from MECT system 10. Image 

reconstructor 34 produces three-dimensional stacks of parallel plane images, or slices, each of 
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which consist of an array of material densities. Some CT images have a resolution of 
512x512x12 bits, and include up to 500 slices in a stack. In a 2D domain, slices can be 
viewed one at a time. An advantage of CT images obtained from image reconstructor 34 is 
that each slice only contains information from one plane. The availability of the stacks of 
parallel data produced by MECT system 10 prompted the development of techniques for 
viewing the volume data as a three-dimensional field, rather than as individual slices. 
Therefore, the volume data can now be viewed from any viewpoint. 

[0041] The method includes repeating 111 step 66 after a patient with metal 
stents and/or metal valves is scanned. Repeating 111 step 66 removes beam-hardening 
artifacts. The beam-hardening artifacts are removed since MECT system 10 scans at multiple 
energies, as described above, to obtain the images of material densities within patient 22 
including the stent. Moreover, repeating 111 step 66 enables looking for restenosis within the 
stent, which is an reappearance of plaque within the stent, by improving the quality of images 
obtained using MECT system 10. It is noted that instead of simulating a scan of the 
phantoms in step 70, an actual placement and an actual scan of the phantoms can be 
performed using MECT system 10. It is also noted that additional lists of materials and 
energies can be created by scanning phantoms with varying compositions of substances, such 
as, for instance, calcium, iodine, and water, at additional energies, such as, third and fourth 
energies. The additional lists can be added to the information in the look-up table. 

[0042] It is noted that the steps of the flowcharts of Figures 3, 4, 5, and 6 can 
be executed in a different order than shown in the Figures. For example, steps 108 and 109 
can be executed before step 106 and after step 104. As another example, steps 106 and 107 
can be performed before step 110 and after step 109. Moreover, it is not necessary to execute 
all procedures identified in Figures 5 and 6 to obtain useful information. 

[0043] Figure 7 illustrates a plurality of images including image 130, image 

132, and image 134. Image 130 is an image reconstructed using conventional CT imaging 

techniques. A high-density organ of image 130 includes calcium 136 and iodine 138 at 

densities that have reconstruct to the same effective CT number or Hounsfeld unit. Hence, 

there is no difference visible between calcium 136 and iodine 138 within image 130. Image 

132 is an iodine density image where iodine 138 is used as a basis material and image 134 is a 

water density image, where water is used as a basis material. Images 132 and 134 were 
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generated by implementing the method of Figures 3, 4, 5, and 6. Specifically, images 132 and 
134 were obtained by implementing step 66, including the basis material decomposition steps 
of the method. The method enables discrimination of calcium 136, visible as a white area, 
and iodine 138, visible as a grey area in image 132. Calcium 136 and iodine 138 can be 
distinguished from each other since each has different attenuation characteristics as a function 
of x-ray spectral energy. Calcium 136 is evident in both images 132 and 134 because it is not 
one of the basis materials used for decomposition. 

[0044] Figure 8 illustrates a plurality of images including image 150, image 
152, and image 154. Image 150 is an image reconstructed using conventional CT imaging 
techniques. High-density organs within image 150 demonstrate beam-hardening artifacts, 
such as non-uniformity, shading and streaking, which occur in some CT images. Image 152 
is an iodine density image where iodine 138 is used as a basis material and image 154 is a 
water density image, where water is used as a basis material. Images 152 and 154 were 
generated by implementing the method of Figures 3, 4, 5, and 6. Specifically, images 152 and 
154 were obtained by implementing the basis material decomposition method in combination 
with step 66. Non-uniformity, shading, and streaking have been removed in images 152 and 
154. The basis material decomposition method in combination with the above described 
implementation of step 66 results in images 152 and 154 which are free of beam hardening 
artifacts. 

[0045] Hence, the herein described systems and methods for detecting 
components of plaque take into account system behavior such as detector response and 
detector signal pileup. The systems and methods have shown improved results as shown in 
Figures 7 and 8. Potential applications of the systems and methods include coronary artery 
imaging and calcium scoring using a protocol specifying contrast injection. Additional 
potential applications include elimination of beam-hardening artifacts in stents and thus 
improvement of image quality within stents, enabling assessment of in-stent restenosis and 
plaque characterization. Moreover, the systems and methods make possible continuous 
screening of high-risk patient populations that may be selected based on genetic or 
environmental history factors. The systems and methods also make possible automatic 
detection of high-risk vascular regions, denoted as unstable plaque regions, based on vascular 
geometry, which is of great benefit to a physician. The systems and methods make possible 
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automatic detection of plaque burden and composition, which is also of benefit to the 
physician. The automatic detection allows for more thorough detection of disease, especially 
early disease markers, in large projection datasets and may permit automated quantification of 
disease at all vascular locations. Automated computation may also be employed to fully 
detect and characterize each plaque region into stages corresponding to different risk groups 
or intervention/therapy groups. The computations may then be highlighted during physician 
review and recommendations for a follow up period and scanning protocols may be made. 
Once cardiovascular disease (CVD) is detected, several methods may be employed to 
automatically detect and diagnose CVD including advanced pattern recognition methods, 
such as, neural networks and model-based algorithms. Progression and regression of CVD is 
also possible if multiple scans are performed by MECT system 10 at varying time points. The 
rate of progression of CVD can be measured in patient 22 and compared to normal 
populations, accounting for age, sex, and other biases. These comparisons will assist in 
identifying patient 22 that is at a high risk of CVD and permit the monitoring of therapy. 
Advanced pattern recognition methods may again be employed on data acquired using MECT 
system 10 and other patient data to ensure that an optimal therapy path is administered for 
patient 22. 

[0046] The herein described systems and methods may also be used as an 
image-based surrogate for drug therapy clinical trials. This allows an assessment of the 
efficacy of pharmacological agents during early drug development and clinical trials. The 
systems and methods can be applied to diagnosis or treatment of a disease, such as CVD, by 
following the regression of the disease. The systems and methods can also be applied to stent 
imaging to assess the presence and extent of re-stenosis. 

[0047] While the invention has been described in terms of various specific 
embodiments which require acquiring projection data at different x-ray spectra, it is possible 
to implement the scheme using only one acquisition if a photon counting MECT is employed 
where the energy of each detected photon is measured. A thresholding method is used in 
which photons having an energy level higher than a threshold are binned as one group and 
photons having an energy level lower than the threshold are binned as another group. In this 
case, an optimal threshold between the upper and lower energy levels is selected from the 
measurements to optimize distinguishing the densities of the basis materials. 
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[0048] While the invention has been described in terms of various specific 
embodiments, those skilled in the art will recognize that the invention can be practiced with 
modification within the spirit and scope of the claims. 
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